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The Bowen- York initial value data typically used in numerical relativity to represent spinning 
black hole are not those of a constant-time slice of the Kerr spacetime. If Bowen- York initial data are 
£ — ' used for each black hole in a collision, the emitted radiation will be partially due to the "relaxation" 

ON , of the individual holes to Kerr form. We compute this radiation by treating the geometry for a 

ON ■ single hole as a perturbation of a Schwarzschild black hole, and by using second order perturbation 

theory. We discuss the extent to which Bowen- York data can be expected accurately to represent 
, Kerr holes. 
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I. INTRODUCTION 



The description of the collision of two black holes, including the total energy radiated and the waveforms to be 



measured by observers far from the collision region, is at this time one of the most active fields of research in general 
relativity. Since the latter theory has a well posed initial value problem, a good deal of effort has been devoted to 
finding solutions for the initial value problem, in the form of initial data sets, that may represent slices of recognizable 
physical processes involving black holes. One of the first examples of this kind was given by Misner JlJ, who derived 
f*** \ an initial data set representing the moment of time symmetry in the head-on collision of two equal mass black holes, 
0^ placed at an arbitrary distance from each other. Further developments along this line have provided a fair number of 
""q I interesting initial data sets that can be interpreted as representing isolated but boosted and/or rotating black 

Q" 1 holes, or collisions involving two or more black holes. 

Once one has the initial data, the next task is to study the evolution and the consequent emission of gravitational 
(^j), waves arising from the collision. Due to the complexity of the evolution equations of general relativity, numerical 
solutions of the full Einstein equations are available at present only for equal mass, nonspinning, holes undergoing a 
. \ head-on (i.e., zero impact parameter) collision. [Q||. 

It has also been observed that for a restricted set of parameters, the evolution can be well described by treating the 
system as a perturbation of a single black hole || . This method, called the "close approximation," has produced results 
that are in remarkable agreement with the full numerical results in the case where the latter are available. An appealing 
feature of the perturbation method is the explicit control over the parameters characterizing the perturbation. This, 
together with the development of some form of "error bars," as in J^,||, can make the method an important tool to 
make predictions of physical processes or to provide comparison cases to test the reliability of numerical methods. 

The initial data set that is usually considered |^,^] in the study of black hole collisions is generated using the confor- 
mal approach, in which one assumes that the spatial metric is conformally flat, and the maximal slicing condition is 
chosen for the extrinsic curvature. In the flat space, it is relatively simple to construct a conformally related extrinsic 
curvature which guarantees that the momentum constraint is solved. The Hamiltonian constraint is subsequently 
solved, either numerically || or via approximations (|J§]. I n this construction there is a series of simplifying assump- 
tions and there is no claim that a "generic" solution has been found. It is therefore not clear that the particular 
solution generated is a true representation of the physical problem in which one is interested. An example of this is 
the Bowen and York (BY) Q solution for a single spinning hole. It is known that this initial data set represents a 
dynamical situation that evolves to a Kerr black hole asymptotically in the future. Initially, however, the spacetime 
is not a Kerr spacetime, but can be thought (somewhat inappropriately) to differ from a Kerr solution in that it has 
some "gravitational wave content," the waves that will be radiated as the spacetime evolves towards Kerr. 

One can argue that this "gravitational wave content" will be radiated in a short time, and that the initial data 
will evolve rapidly to a stationary Kerr black hole configuration, and therefore will not greatly affect the radiation 
produced in a black hole collision, as long as the holes are released far from each other. The question is important 
enough to deserve a more careful answer. In addition, numerical relativity codes cannot be accurately run for long 
evolution times, so initial data will have to be specified at fairly late times, with holes fairly close together, and with 
the possibility that the "gravitational wave content" of the initial data is a major part of the outgoing radiation. 
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The purpose of this paper is to analyze the radiation from the individual holes. More specifically, we use the theory 
of perturbations of the Schwarzschild spacetime. We consider that we have a family of spacetimes depending on the 
parameter e, and that, in appropriate coordinate systems, metrics of the family can be expanded as 

9»u=9t) + eh$+e*h%+... . (1) 

Here gjS) the Schwarzschild metric, fify is called the first order perturbation, and hffi is called the second order 
perturbation. To analyze the "Bowen and York spacetime," the spacetime that evolves from Bowen and York initial 
data for a single spinning hole, we choose as the expansion parameter the angular momentum J, and we analyze both 
the Kerr solution and the BY spacetime to second order in J. The second order expansions are then compared and 
we find that we need only to evolve the difference between the BY spacetime and the Kerr spacetime. 

The organization of this paper is as follows. In section II we give a perturbation analysis of the BY initial data. 
The method of evolving this initial data is described in section III. In section IV the results for radiation emitted are 
presented and discussed. In an Appendix we show how the Kerr metric can be written as an expansion, in angular 
momentum, about the Schwarzschild spacetime. 

When it is useful to specify orders of expansion, and/or multipole indices, we shall use a leading subscript to denote 
the multipole index £, and a superscript, in parenthesis, to the right of a perturbation variable, to indicate order in 
J. The quantity oip^ 2 \ for example, is a monopole perturbation, second order in J. 



II. THE BOWEN-YORK SINGLE ROTATING BLACK HOLE 

The Bowen- York Q construction of initial data assumes that space-time contains a (constant t) slice, where the 
3- metric can be written in the form — ip 4 fij, where is the metric for the flat background ds 2 = dR 2 + R 2 (dd 2 + 
sin 2 6 dcj) 2 ), while the extrinsic curvature is given by K%j = ij)~ 2 Kij, and satisfies K\ — 0. With these assumptions the 
initial value constraint equations take the form 



and 



K v l3 = 0, (2) 



v 2 v = -\k t] k i ^-\ (3) 



where the Laplacian, and the covariant differentiation (denoted by a vertical bar) are taken with respect to the flat 
space metric. The problem is not fully specified until appropriate boundary conditions are imposed on ip. The 
Bowen- York prescription for a single black hole is that, given a certain constant a, Eq. (|^) holds for R > a, and 

H + ^ = 0forE = a , (4) 



and 



ip > , lim V = 1- (5) 

R — >oo 



The particular solution we are interested in corresponds to Q 



K = [(JAR)®R + R® (JAR)], (6) 

where R is the "position" vector in the flat space background, and J is a vector constant. Without loss of generality, 

we choose J = Jk, where k is a unit vector pointing along (the positive direction of) the polar axis of the R, 9, <f> 
coordinate system. With this choice, the only nonvanishing components of K are 

K R< p = K^r = ^ sin 2 9 , (7) 

We then find 

i^ir j = 18^sin 2 0. (8) 
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We now need to solve Eq. (^|). For general values of </, this can only be done numerically. However, we are really 
interested in solutions near J = ("slow rotation") and an expansion in powers of J is appropriate. The zeroth order 
equation is 

V 2 V = , (9) 

and the solution that satisfies the boundary conditions is 

V> (0) = 1 + a/R . (10) 

The lowest order correction to tp may now be constructed by linearizing tfj about in Eq. (^). If we formally 
write 

i> = v^ (0) + J 2 ^ {2) + ■■■ , (ii) 

the resulting equation is 

V V > = -ji, sin* «(l + (12 ) 

and the sin 2 9 factor may be expanded in Legendre polynomials, as 2[Po(cos6*) — P2(cos 6>)]/3. 
We next write i[)( 2 '(R, 9) as 

Tp (2 HR,9) = o^ 2) (R)Po(co8 9) + 2 ^ {2 \R) P 2 (cos9) (13) 
and we find that oV'^H-^) an d 2ip^ 2 \R) satisfy the equations 

rf 2 [o^)] | 2 d[of'] _ 3 

and 



dR 2 ' R dR 2i?6 V 1 + Rj ' 



ft^l 2£t^i a 2 ^)_J_r i + «r 7 d5) 



The solutions of these equations, satisfying the boundary conditions are 



and 



(2) a 4 + R 4 + 5aR(R + a) 2 
^ = iOa^R + af (16) 



^ i2) = -^Srw- (17) 

With Eq. ([l6]) we see that the conformal factor has the form 

* = l + ± + ^ + 0(1/1?) + 0(J>). (18) 

But by the definition of the ADM mass (see for instance f|), i> = 1 + Madm/(2-R) + 0(l/i? 2 ), and hence the ADM 
mass, to second order in J, is 

J2 

Madm = 2a+ — — . (19) 
20a c> 

This dependence of Madm on J, for constant a, is plotted in Fig.|l| as the dashed curve. Two other computations of 
-Madm are also plotted, based on the integral P,|l0||, 



Madm = ^- / — y^dv + - / ^smOdO . (20) 
32tt J r>a V 7 2 J 
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FIG. 1. Comparison of approximate solutions to the initial value problem with a full numerical integration performed with 
a multigrid method. Our code uses a square grid in compactified coordinates p = R/(a + R), 6 with 200 x 200 grid points. It 
achieves an accuracy similar to the test runs of Choptuik and Unruh (lo| . The figure at the left compares the ADM mass of 
approximate and numerical solutions. The two approximations are described in the text. The figure on the right shows the 
percentile difference between the full numerical solution for tp and the second order approximation, for J/a 2 — 4, as a function 
of p, 8. We see that the approximation works very well. 



The solid curve shows the result of substituting tp, to second order in J, into the integrals in Eq. (po|). The result 
agrees with Eq. (|l9) to second order in J, but is significantly smaller for large J. (The difference is due to terms 
higher order in J; the expressions are identical up to second order). The dark points show the result of a multigrid 
numerical code we wrote to solve the fully nonlinear initial value problem of general relativity with axisymmetry, 
similar to that used by Choptuik and Unruh jul . 

The results presented to this point are formally expansions to second order in the dimensionless parameter J/a 2 , 
but a itself is not a physical parameter, so the physical meaning of this expansion is unclear. For that reason we 
now convert our result to an expansion in the parameter J /M 2 . (Here, and in subsequent expressions, we drop the 
"ADM" subscript on M, The symbol M will always represent the ADM mass.) We then consider in the expressions 
we gave for the conformal factor, 



a 
R 



J 2 (o^+ 2^2 P 2 (COS 6)) 



(21) 



that the parameter a is really given by a{M, J). One can obtain explicit formulas for this by considering the expression 
for a(M, J) obtained by inverting the approximate expression Eq. ([l9]). Alternatively, one can keep the implicit 
dependence of a on J, and can replace it in the last step of a calculation, using tabulated values for a(M, J) obtained 
from the multigrid numerical code. In this case the "second order part" is the part of ip that remains after the part 
zeroth order in J/M 2 is subtracted. 

If the relation of a and J is taken from Eq. (|l9|); the explicit formulas (correct to second order) are, 



i> = *(°> + (J/M 2 ) 2 ^ = *(°> + (J/M 2 ) 2 ( * (2) + 2 * (2) P 2 (cos#)) 



with 



#(°) = 1 + M/2R 



(22) 



(23) 



M 1 + (5M/2i?))(l + M/2R) 2 + (M/2R) i - (1 + M/2Rf 



R 



5(1 + M/2R) 5 



(M/Rf 
10(1 + M/2Rf 



(24) 
(25) 
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To obtain the background metric in the usual Schwarzschild form we introduce the radial coordinate r with 

R= (r- M + \jr 2 - 2Mr) /2 . (26) 
Since the conformal factor, up to second order, can be written as 

j \ 2 



v> 4 = (V°>) 



1 + 4 



M 2 I *(°) 



the 3-geometry to second order is given by 



2A/ N 
1 )g r 

r 



r 2 r 2 sin 2 I 



and, in terms of the r variable, we may write 



M 3 



\f(o) lOr 3 



\M 2 ) *(0) 



(27) 



(28) 



(29) 



III. PERTURB ATIVE EVOLUTION OF THE BOWEN-YORK INITIAL DATA 



We adopt the notation of Regge and Wheeler (TTJ for the separation into parities and the multipole decomposition 
of the perturbations. Because we are considering axisymmetric situations, all multipole decompositions are given in 
terms of Legendre polynomials Pi(cos6). On the initial constant-i hypersurface of the BY spacetime, we can read off 
the perturbations of the 3-geometry from Eq. ( p8| ) . The perturbations are purely second order even parity, and in the 
Regge- Wheeler notation, are: 



BY 



K 



BY 



A(J/M 2 ) 2 2 *( 2 )/*(°) = -4J 2 /(10Mr 3 ) 



G 



BY 



BY 



. 



(30) 



The first order perturbations contained in the BY spacetime are those specified by the extrinsic curvature in Eq. 
(0), which can be reexpressed as 



dR 



K rcl , = i/j K r< j, =ip -7- K R<f> = ~, p. kttt- 

1 dr I v r 2 yjl - 2M/r 



3J 



sin 



(31) 



and is identical to the extrinsic curvature given in Eq. ( pq ) for the Kerr geoemtry. The explicit form of the metric 
perturbations, to first order in J, depends on the gauge (first order coordinate fixing) we choose. Let us choose the 
coordinates to first order so that the initial BY metric is the same as the first order Kerr metric given by Eq. fl5^) of 
the Appendix. That is, let us choose 



(32) 



to be the only nonvanishing first order initial perturbation. Since this perturbation is purely I = 1, the Einstein 
equations require that there be no gauge independent time variation in this perturbation. Let us choose, therefore, 
to have the perturbation in Eq. (|32| ) be the first order perturbation for all time. This is equivalent to the statement 
that we are choosing coordinates, to first order in J, so that the BY and the Kerr spacetimes agree, at all times, to 
first order in J. 

In principle, the radiation in the BY spacetime would be found from the quadrupolar part of Einstein's equations 
to second order in J/M 2 . These equations contain terms linear in the 1 = 2 second order metric perturbations, and 
quadratic in the first order perturbations. The general structure of these equations is discussed in Ref. @. In that 



reference it is shown how these second order equations can be combined into an equation like that of Zerilli J_2 13 
for first order equations. The second order equivalent of the Zerilli equation differs only in that it has "source" 
terms quadratic in the first order perturbations. Since we know the first order perturbations, for all time, for the BY 
spacetime, we know the source term. The wave equation can therefore be solved numerically and the radiation signal 
found from the solution for the wave variable at large radius. 



5 



In practice, this computation can be made much simpler. Since the first order perturbations in the BY and the 
Kerr spacetimes are identical, the source terms will be identical in the second order Einstein equations for BY and 
for Kerr. We can exploit this, by evolving only the difference between the second order BY and Kerr perturbations, 
following the technique used by Cunningham et al. jL4|. To do this, we take as our second order £ = 2 variable the 
Moncrief |H| wave variable 

(2r-4rM) ( 2 d[ 2 K^] 2 d[ 2 G^} (2) 19 ,(2)\ . „ (a) r _„. 

2XM(r,t) = — - I -2r — — — +6r r + 2r 2 H 2 > -12 2 h\ j + ir 2 K^ > (33) 

in which 2 K^ 2 \ 2 G^ 2 \ 2 h^ and 2 H 2 2 ^ are second order £ = 2 perturbations. This Moncrief wave variable has two 
very useful features. First, 2 \m is constructed only from perturbations in the 3-geometry on the initial hypersur- 
faces. Second, it is invariant under second order coordinate transformations, that is, under transformations of form 
(x») New = (x^) old + in which f is second order. 

We are using here the same normalization for our wave function as in Ref. Q. This normalization is formally the 
same as that of Zerilli |l^ ], except that we expand in Pi rather than Yf m ; as a result our 2 xm is related to the variable 
K 2 m of Zerilli |l2| by 2 \m = ^/5/47ri^ 2 M- For a discussion of the relationship of the Moncrief and Zerilli variables, 
and various normalizations, see Ref. It is straightforward to generalize the analysis to arbitrary I > 2. 

We use the Moncrief wave functions 2 Xyt to describe the Bowen York spacetime, and 2XM° rr f° r Kerr. The initial 
value for 2 Xm Y are taken from Eq. j30| ) and the value of 2 XM° rr , initially and for all time, are taken from the expansion 
in the Appendix, from which we have, 

4J2 

2jff ° KC1T = "3A/r2(r-2M)' (34) 
2 Hf err = (35) 



2J 2 
3MV 2 



^ = 3W (36 ^ 
Kcrr 2J 2 (4M + r) 

2K = (3?) 

Korr _ J\2M + r) 

2 ° - 3MV ■ (38) 

We now define a "radiative" Moncrief wave function by 

, rad — ,,By _Kerr fon\ 

Since the first order perturbations, and therefore the source terms, for 2 Xm^ and for 2Xm git are ^ ne same, 2X rad 
satisfies a homogeneous Zerilli function 

d 2 [ 2X rad (t,r)} d 2 [ 2X ™ d (t,r)} rad 



dt 2 



y(r*) 2 X rad (i,O=0, (40) 



where 



and where V is the t = 2 Zerilli potential 



2Mln[r/(2M) - 1] , (41) 



V(r) =6(1-2^ 4r3+4r ^ + 6 ;5 2 +3M3 ■ (42) 



r J r 3 (2r + 3M) 2 

A word of explanation is appropriate about the properties of 2 xm under a coordinate transformation. As already 
stated, 2 xm is invariant under a transformation in which the coordinates change only to second order. This is why 
we can use Eqs. (|i|) - ( p8[ ) to evaluate 2 Xm", and Eq. ( |30| ) for 2 Xm Y ' though the second order metric perturbations 
used are clearly in different gauges. The point is that we know that the second order perturbations can be brought 
into the same second order gauge with second order gauge transformations, and that this has no effect on 2 x^ crr or 
on 2Xm F i and hence no effect on 2 x rad - It should be realized that 2 xm is not invariant under a first order change of 
coordinates. If we were to perform, say, a first order change in the coordinates used in the Appendix, then the value 
of 2XM 6rr we would compute would change. It is important, therefore, that no first order change in coordinates is 
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needed for the Kerr expansion in the Appendix, or the BY spacetime in Sec.[n]. They are already in the same first 
order gauge. 

The initial 2XM for this equation is simply the known difference between the initial forms of 2Xm^ ana - 2XM° rr ' ana - 
turns out to be 



, 2 J z ( -5r d - 7 Mr* + 2hM A r + 60M d 



l*=o = 5jy^3 [ 27^3M ) ■ (43) 

There are no second order perturbations to the extrinsic curvature of a constant time slice of Kerr, or in the BY initial 
data. The first of these conclusions follows from an explicit computation based on the metric in the Appendix. One 
finds that the extrinsic curvature contains only odd powers of J. (This conforms to the intuition that suggests that 
reversing the direction of J should reverse the sign of extrinsic curvature.) The conformally related extrinsic curvature 
Kij for the BY initial data is given to all orders by Eq. (||). The second order perturbations in the conformal factor 
mean that the extrinsic curvature = ip~ 2 Kij will contain perturbations of odd orders in J, due to the perturbations 
of even order contained in tp. Since there are no second order contributions to the extrinsic curvature, of cither BY 
or Kerr, it follows that 



|(rf-) 



= . (44) 

t=o 



The wave equation of Eqs. (|4(]) - (|4^), with the cauchy data of Eqs. (f43| ) - (Q), is simply solved numerically for 
2X rad (^ r )j an d from the solution we can find 

2Xm \t,r) = 2X rad (t,r) + 2X ^ I (r) , (45) 

where 2XM err (^) is the known, time independent, Kerr solution. The radiative power (see Ref. |Q) contained in the 
BY spacetime is then given by 

p — Hi^y-Hi^y ■ <46) 

and this is the gravitational radiation power emitted as the BY solution settles into its Kerr final form. The energy 
radiated is the time integral of this expression. 
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FIG. 2. The radiated waveform for a single spinning Bowen and York black hole, computed treating the spacetime as a 
perturbation of Schwarzschild. The waveform, extracted at r — 20M, is dominated by quasinormal ringing which persists long 
enough to contaminate radiation coming out from a head on collision of black holes. We display the waveforms corresponding to 
the different choices of background mass used in perturbation theory. The significant differences among these curves illustrates 
the importance of the choice of the background mass in perturbative calculations, i.e., on the way in which the ADM mass M 
is taken to depend on J, for a given value of a. The bold curve uses the "correct" ADM mass computed numerically with the 
multigrid code. The thin curve was computed with expression M = 2a + J 2 /(20a 3 ) (see text) which is correct to second order 
in J. The dashed curve uses M = 2a, and ignores the influence of J on M, for given a. Previous experience with perturbation 
calculations for which numerical relativity comparisons were available strongly suggests that the curve with the numerically 
computed ADM mass is the most accurate. 



IV. RESULTS AND DISCUSSION 

In Fig. ^| we show the radiated waveform, as a function oft, for fixed, large r, from which we may infer the effective 
time for the decay of the Bowen- York rotating black hole into its final Kerr state. It is clear that the wave form is 
dominated by quasinormal ringing. This means that the "initial burst" of energy generated as the hole relaxes to 
Kerr form can contaminate the evolution for some time. 

In Fig. U we show the total energy radiated as a function of J/M 2 . Whatever choice we make for the J dependence 
of the ADM mass, our perturbation calculation for x rad is formally correct only to second order in J/M 2 . Since the 
radiated energy is quadratic in x rad , the energy results displayed in Fig.^ are formally correct only to fourth order in 
J/M 2 , the lowest nontrivial order. The results in Fig.||, cannot therefore be trusted for J/M 2 near the astrophysically 
interesting limit J/M 2 « 1. We suspect that the curve corresponding to the numerical ADM mass is reliable within 
a factor of two or so up to J/M 2 — 0.8. A more accurate evaluation will require either fully nonlinear numerical 
relativity, or a calculation using second order perturbations around the Kerr solution. 

According to Fig.||, the "BY relaxation energy," the energy emitted as a single BY hole relaxes to a Kerr hole 
appears to be small. It should be kept in mind, however, that the total radiation in a black hole coalescence can 
be comparably small. For a head on collision of nonspinning holes the total radiated energy is of order 10~ 3 of the 
total ADM mass. Head on collisions, of course, are not of primary astrophysical interest. For the "merger" phase of 
equal mass holes, radiated energy is expected to be several percent JL7| . In this case, the BY relaxation energy would 
be negligibly small. It would, furthermore, be emitted within a few quasinormal periods of a single hole, while the 
merger and ringdown of the final hole formed would require a time an order of magnitude longer. 
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FIG. 3. The total radiated energy computed for "BY relaxation." The method is formally correct only to lowest (fourth) 
order in the expansion parameter J/M\ DM . For a description of the various ADM masses involved, see the caption of figure 2. 



the Pennsylvania State University and its Office for Minority Faculty Development, and the Eberly Family Research 
Fund at Penn State. We also acknowledge support of CONICET and CONICOR (Argentina). JP also acknowledges 
support from the Alfred P. Sloan Foundation through an Alfred P. Sloan fellowship. 



[1] C. Misner, Phys. Rev. D118, 1110 (1960). 

[2] J. Bowen, J. York, Phys. Rev. D21, 2047 (1980). 

[3] G. B. Cook, M. W. Choptuik, M. R. Dubai , Phys. Rev. D 47, 1471 (1993); G. B. Cook, Phys. Rev. D 50, 5025 (1994); 

Ph.D. thesis, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, 1990. 
[4] P. Anninos, D. Hobill, E. Seidel, L. Smarr, W.-M. Suen, Phys. Rev. Lett. 71, 2851 (1993); Phys. Rev. D 52, 2044 (1995). 
[5] J. Baker, A. Abrahams, P. Anninos, S. Brandt, R. Price, J. Pullin, E. Seidel, Phys. Rev. D55, 829 (1997). 
[6] R. Price, J. Pullin, Phys. Rev. Lett. 72, 3297 (1994). 

[7] R. Gleiser, O. Nicasio, R. Price, J. Pullin, Class. Quan. Grav. 13, L117 (1996). 

[8] R. Gleiser, O. Nicasio, R. Price, J. Pullin, Phys. Rev. Lett. 77, 4483 (1996). 

[9] J. Pullin, Fields Inst. Commun. 15, 117 (1997). 
[10] M. Choptuik, W. Unruh, Gen. Rel. Grav. 18, 813 (1986). 
[11] T. Regge, J. Wheeler, Phys. Rev. 108, 1063 (1957). 
[12] F. J. Zerilli, Phys. Rev. Lett. 24 737 (1970). 
[13] F. J. Zerilli, Phys. Rev. D2 2141 (1970). 

[14] C. Cunningham, R. Price, V. Moncrief, Astroph. J. 236, 674 (1980). 

[15] V. Moncrief, Ann. Phys. (NY) 88, 323 (1974). 

[16] CO. Lousto and R.H. Price, Phys. Rev. D55, 2124 (1997). 

[17] E.E. Flanagan and S.A. Hughes, preprint gr-qc 9701039. 



9 



APPENDIX: THE KERR METRIC AS A PERTURBATION OF A SCHWARZSCHILD BLACK HOLE 



The Kerr metric in Boyer - Lindquist coordinates takes the form 



ds 2 = 



A - ( JIM) 2 sin 2 6 



dt 2 - 2( J/M) sin 2 9 



p 2 + ( J/M) 2 - A 



{p 2 + (J/M) 2 ) 2 - A (J/M) 2 sin 2 9 2 



A 



dp 2 



dtd<j> 
£d6> 2 



(47) 



where £ = p 2 + (J/M) 2 cos 2 and, A = p 2 + (J/M) 2 - 2Mp. 

If we assume J < M 2 , @ is defined only for p > M + ^M 2 - (J/M) 2 , since A = for p = M + ^M 2 - (J/M) 2 . 

The metric ( fi"7| ) reduces to the Schwarzschild metric, in the range p > 2M, for J = 0. It seems reasonable, therefore, 
to try to find an expansion of (^), in powers of J, as a perturbation of a Schwarzschild black hole. This expansion, 
however, would fail near p = 2M, because the metric coefficient g pp does not have the required analyticity properties. 
To avoid this problem we introduce a new coordinate r, such that 



r 2 - 2Mr = p 2 + (J/M) 2 - 2Mp 



(48) 



With this definition we have that A = corresponds to r = 2M. We may invert (48) to 



p = M + y/(r - M) 2 - (J/M) 2 (49) 
Then, for r > 2M, and J < M 2 , the right hand side may be expanded in powers of J. The leading terms are 

J 2 



p = r 



2M 2 (r-M) 



(50) 



In fact, one can easily check that all the metric coefficients admit a convergent power series expansion in J. To leading 
order we have 



9tt 



= 1- 



2M 



-1 + 



J 2 (r + 2M) 



r 

7 ^-1 

2- sin0— Pi 

r o0 



2M 
r 



3Mr 2 (r-M)(r-2M) 
0QJ/M 2 ] 3 ) 

J 2 (Mr + r 2 +M 2 ) 



4J 2 



3Mr 2 (r - 2M) 



P 2 + 0([J/M 2 



3M 2 r 2 (r - M) 2 

J 2 (2M 2 + r 2 ) 2J 2 (AM- 
Po - 



-Po 



2J 2 



3M 2 r 2 



P 2 + 0([J/M 2 



r 2 sin 2 i 



3M 2 r 3 (r - M) 

J 2 (2M 2 + r 2 ) 



Po 



3M 2 r 3 
2J 2 (4M 



r) J 2 (r + 2M) d 2 

— Pi rTTil ^ -7^P2 



3M 2 r 3 d9 2 ' 



0([J/M 2 f) 



J 2 (r + 2M) d 
SM 2 ^ ~d9 



3M 2 r 3 (r-M) 3M 2 r 3 
0([J/M 2 } 4 ) 



-Pi 



(51) 
(52) 
(53) 
(54) 

(55) 



where P = 1, Pi = cos# and, P2 = (3/2) cos 2 9 — 1/2 are Legendre polynomials. 

From this metric it is straightforward to compute, to first order in J, the extrinsic curvature of a t — constant 
surface. If we let n be the future directed normal to a t = constant hypersurface, then, the extrinsic curvature is 
Kij — — Jiju, where the bar denotes covariant differentiation with respect to the 3-geometry. The normal has only 

a single covariant component n t which, to first order in J, is rit = — — 2M Jr. With this, a straightforward 
computation shows that the only nonvanishing first order components of are 



K r 



K, 



3J 



4>r 



rVl - 2M/r 



sin 2 6 



(56) 
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